Shaping the future of preclinical development of successful disease-modifying drugs against Alzheimer's disease: a systematic review of tau propagation models.

The transcellular propagation of the aberrantly modified protein tau along the functional brain network is a key hallmark of Alzheimer's disease and related tauopathies. Inoculation-based tau propagation models can recapitulate the stereotypical spread of tau and reproduce various types of tau inclusions linked to specific tauopathy, albeit with varying degrees of fidelity. With this systematic review, we underscore the significance of judicious selection and meticulous functional, biochemical, and biophysical characterization of various tau inocula. Furthermore, we highlight the necessity of choosing suitable animal models and inoculation sites, along with the critical need for validation of fibrillary pathology using confirmatory staining, to accurately recapitulate disease-specific inclusions. As a practical guide, we put forth a framework for establishing a benchmark of inoculation-based tau propagation models that holds promise for use in preclinical testing of disease-modifying drugs.

Hiding in plain sight: Complex interaction patterns between Tau and 14-3-3ζ protein variants.

Tau protein is an intrinsically disordered protein that plays a key role in Alzheimer's disease (AD). In brains of AD patients, Tau occurs abnormally phosphorylated and aggregated in neurofibrillary tangles (NFTs). Together with Tau, 14-3-3 proteins - abundant cytosolic dimeric proteins - were found colocalized in the NFTs. However, so far, the molecular mechanism of the process leading to pathological changes in Tau structure as well as the direct involvement of 14-3-3 proteins are not well understood. Here, we aimed to reveal the effects of phosphorylation by protein kinase A (PKA) on Tau structural preferences and provide better insight into the interaction between Tau and 14-3-3 proteins. We also addressed the impact of monomerization-inducing phosphorylation of 14-3-3 at S58 on the binding to Tau protein. Using multidimensional nuclear magnetic resonance spectroscopy (NMR), chemical cross-linking analyzed by mass spectrometry (MS) and PAGE, we unveiled differences in their binding affinity, stoichiometry, and interfaces with single-residue resolution. We revealed that the interaction between 14-3-3 and Tau proteins is mediated not only via the 14-3-3 amphipathic binding grooves, but also via less specific interactions with 14-3-3 protein surface and, in the case of monomeric 14-3-3, also partially via the exposed dimeric interface. In addition, the hyperphosphorylation of Tau changes its affinity to 14-3-3 proteins. In conclusion, we propose quite complex interaction mode between the Tau and 14-3-3 proteins.

Conformation and Affinity Modulations by Multiple Phosphorylation Occurring in the BIN1 SH3 Domain Binding Site of the Tau Protein Proline-Rich Region.

An increase in phosphorylation of the Tau protein is associated with Alzheimer's disease (AD) progression through unclear molecular mechanisms. In general, phosphorylation modifies the interaction of intrinsically disordered proteins, such as Tau, with other proteins; however, elucidating the structural basis of this regulation mechanism remains challenging. The bridging integrator-1 gene is an AD genetic determinant whose gene product, BIN1, directly interacts with Tau. The proline-rich motif recognized within a Tau(210-240) peptide by the SH3 domain of BIN1 (BIN1 SH3) is defined as 216PTPP219, and this interaction is modulated by phosphorylation. Phosphorylation of T217 within the Tau(210-240) peptide led to a 6-fold reduction in the affinity, while single phosphorylation at either T212, T231, or S235 had no effect on the interaction. Nonetheless, combined phosphorylation of T231 and S235 led to a 3-fold reduction in the affinity, although these phosphorylations are not within the BIN1 SH3-bound region of the Tau peptide. Using nuclear magnetic resonance (NMR) spectroscopy, these phosphorylations were shown to affect the local secondary structure and dynamics of the Tau(210-240) peptide. Models of the (un)phosphorylated peptides were obtained from molecular dynamics (MD) simulation validated by experimental data and showed compaction of the phosphorylated peptide due to increased salt bridge formation. This dynamic folding might indirectly impact the BIN1 SH3 binding by a decreased accessibility of the binding site. Regulation of the binding might thus not only be due to local electrostatic or steric effects from phosphorylation but also to the modification of the conformational properties of Tau.

Biomolecular Complexation on the "Wrong Side": A Case Study of the Influence of Salts and Sugars on the Interactions between Bovine Serum Albumin and Sodium Polystyrene Sulfonate.

In the protein purification, drug delivery, food industry, and biotechnological applications involving protein-polyelectrolyte complexation, proper selection of co-solutes and solution conditions plays a crucial role. The onset of (bio)macromolecular complexation occurs even on the so-called "wrong side" of the protein isoionic point where both the protein and the polyelectrolyte are net like-charged. To gain mechanistic insights into the modulatory role of salts (NaCl, NaBr, and NaI) and sugars (sucrose and sucralose) in protein-polyelectrolyte complexation under such conditions, interaction between bovine serum albumin (BSA) and sodium polystyrene sulfonate (NaPSS) at pH = 8.0 was studied by a combination of isothermal titration calorimetry, fluorescence spectroscopy, circular dichroism, and thermodynamic modeling. The BSA-NaPSS complexation proceeds by two binding processes (first, formation of intrapolymer complexes and then formation of interpolymer complexes), both driven by favorable electrostatic interactions between the negatively charged sulfonic groups (-SO3-) of NaPSS and positively charged patches on the BSA surface. Two such positive patches were identified, each responsible for one of the two binding processes. The presence of salts screened both short-range attractive and long-range repulsive electrostatic interactions between both macromolecules, resulting in a nonmonotonic dependence of the binding affinity on the total ionic strength for both binding processes. In addition, distinct anion-specific effects were observed (NaCl < NaBr < NaI). The effect of sugars was less pronounced: sucrose had no effect on the complexation, but its chlorinated analogue, sucralose, promoted it slightly due to the screening of long-range repulsive electrostatic interactions between BSA and NaPSS. Although short-range non-electrostatic interactions are frequently mentioned in the literature in relation to BSA or NaPSS, we found that the main driving force of complexation on the "wrong side" are electrostatic interactions.

A new fibrillization mechanism of ß-lactoglobulin in glycine solutions.

Even though amyloid aggregates were discovered many years ago the mechanism of their formation is still a mystery. Because of their connection to many of untreatable neurodegenerative diseases the motivation for finding a common aggregation path is high. We report a new high heat induced fibrillization path of a model protein ß-lactoglobulin (BLG) when incubated in glycine instead of water at pH 2. By combining atomic force microscopy (AFM), transmission emission microscopy (TEM), dynamic light scattering (DLS) and circular dichroism (CD) we predict that the basic building blocks of fibrils made in glycine are not peptides, but rather spheroid oligomers of different height that form by stacking of ring-like structures. Spheroid oligomers linearly align to form fibrils by opening up and combining. We suspect that glycine acts as an hydrolysation inhibitor which consequently promotes a different fibrillization path. By combining the known data on fibrillization in water with our experimental conclusions we come up with a new fibrillization scheme for BLG. We show that by changing the fibrillization conditions just by small changes in buffer composition can dramatically change the aggregation pathway and the effect of buffer shouldn't be neglected. Fibrils seen in our study are also gaining more and more attention because of their pore-like structure and a possible cytotoxic mechanism by forming pernicious ion-channels. By preparing them in a simple model system as BLG we opened a new way to study their formation.

Phosphorylated and Phosphomimicking Variants May Differ-A Case Study of 14-3-3 Protein.

Protein phosphorylation is a critical mechanism that biology uses to govern cellular processes. To study the impact of phosphorylation on protein properties, a fully and specifically phosphorylated sample is required although not always achievable. Commonly, this issue is overcome by installing phosphomimicking mutations at the desired site of phosphorylation. 14-3-3 proteins are regulatory protein hubs that interact with hundreds of phosphorylated proteins and modulate their structure and activity. 14-3-3 protein function relies on its dimeric nature, which is controlled by Ser58 phosphorylation. However, incomplete Ser58 phosphorylation has obstructed the detailed study of its effect so far. In the present study, we describe the full and specific phosphorylation of 14-3-3ζ protein at Ser58 and we compare its characteristics with phosphomimicking mutants that have been used in the past (S58E/D). Our results show that in case of the 14-3-3 proteins, phosphomimicking mutations are not a sufficient replacement for phosphorylation. At physiological concentrations of 14-3-3ζ protein, the dimer-monomer equilibrium of phosphorylated protein is much more shifted towards monomers than that of the phosphomimicking mutants. The oligomeric state also influences protein properties such as thermodynamic stability and hydrophobicity. Moreover, phosphorylation changes the localization of 14-3-3ζ in HeLa and U251 human cancer cells. In summary, our study highlights that phosphomimicking mutations may not faithfully represent the effects of phosphorylation on the protein structure and function and that their use should be justified by comparing to the genuinely phosphorylated counterpart.

Quantitation of Human 14-3-3ζ Dimerization and the Effect of Phosphorylation on Dimer-monomer Equilibria.

14-3-3 proteins are universal regulatory proteins and their function depends on their oligomeric form which may alter between the monomeric, homodimeric and heterodimeric states. The populations of individual oligomeric forms are controlled by Kd values of the dimer-monomer equilibria between the involved isoforms. This complex picture is extended by post-translational modifications, e.g. phosphorylation. In this work, we describe the equilibria between monomers, homo- and heterodimers of the 14-3-3ζ isoform in the unmodified and phosphorylated form. To cover a wide range of dimerization affinities, we combined solution NMR, microscale thermophoresis, native PAGE, and a set of novel fluorescence assays. Using a FRET based assay, we also determined the kinetic parameters of dimerization. We found that phosphorylation of 14-3-3ζ at Ser58 increases its homodimeric Kd value by 6 orders of magnitude. The presented assays allow to efficiently monitor 14-3-3ζ dimerization as a function of external factors, such as temperature, salt concentration, and client protein binding. For instance, we obtained values of both transient and equilibrium thermodynamic constants for the dimerization, and observed a substantial decrease of 14-3-3ζ dimer dissociation rate upon binding to the doubly phosphorylated regulatory domain of tyrosine hydroxylase. In summary, our work provides a conceptual framework to characterise the isoform exchanges of homo- and heterodimers which can significantly deepen our knowledge about the regulatory function of 14-3-3 proteins.

NMR Studies of Tau Protein in Tauopathies.

Tauopathies, including Alzheimer's disease (AD), are the most troublesome of all age-related chronic conditions, as there are no well-established disease-modifying therapies for their prevention and treatment. Spatio-temporal distribution of tau protein pathology correlates with cognitive decline and severity of the disease, therefore, tau protein has become an appealing target for therapy. Current knowledge of the pathological effects and significance of specific species in the tau aggregation pathway is incomplete although more and more structural and mechanistic insights are being gained using biophysical techniques. Here, we review the application of NMR to structural studies of various tau forms that appear in its aggregation process, focusing on results obtained from solid-state NMR. Furthermore, we discuss implications from these studies and their prospective contribution to the development of new tauopathy therapies.

Universal two-point interaction of mediator KIX with 9aaTAD activation domains.

The nine-amino-acid activation domain (9aaTAD) is defined by a short amino acid pattern including two hydrophobic regions (positions p3-4 and p6-7). The KIX domain of mediator transcription CBP interacts with the 9aaTAD domains of transcription factors MLL, E2A, NF-kB, and p53. In this study, we analyzed the 9aaTADs-KIX interactions by nuclear magnetic resonance. The positions of three KIX helixes α1-α2-α3 are influenced by sterically-associated hydrophobic I611, L628, and I660 residues that are exposed to solvent. The positions of two rigid KIX helixes α1 and α2 generate conditions for structural folding in the flexible KIX-L12-G2 regions localized between them. The three KIX I611, L628, and I660 residues interact with two 9aaTAD hydrophobic residues in positions p3 and p4 and together build a hydrophobic core of five residues (5R). Numerous residues in 9aaTAD position p3 and p4 could provide this interaction. Following binding of the 9aaTAD to KIX, the hydrophobic I611, L628, and I660 residues are no longer exposed to solvent and their position changes inside the hydrophobic core together with position of KIX α1-α2-α3 helixes. The new positions of the KIX helixes α1 and α2 allow the KIX-L12-G2 enhanced formation. The second hydrophobic region of the 9aaTAD (positions p6 and p7) provides strong binding with the KIX-L12-G2 region. Similarly, multiple residues in 9aaTAD position p6 and p7 could provide this interaction. In conclusion, both 9aaTAD regions p3, p4 and p6, p7 provide co-operative and highly universal binding to mediator KIX. The hydrophobic core 5R formation allows new positions of the rigid KIX α-helixes and enables the enhanced formation of the KIX-L12-G2 region. This contributes to free energy and is the key for the KIX-9aaTAD binding. Therefore, the 9aaTAD-KIX interactions do not operate under the rigid key-and-lock mechanism what explains the 9aaTAD natural variability.

Cysteine Modification by Ebselen Reduces the Stability and Cellular Levels of 14-3-3 Proteins.

The 14-3-3 proteins constitute a family of adaptor proteins with many binding partners and biological functions, and they are considered promising drug targets in cancer and neuropsychiatry. By screening 1280 small-molecule drugs using differential scanning fluorimetry (DSF), we found 15 compounds that decreased the thermal stability of 14-3-3ζ Among these compounds, ebselen was identified as a covalent, destabilizing ligand of 14-3-3 isoforms ζ, ε, γ, and η Ebselen bonding decreased 14-3-3ζ binding to its partner Ser19-phosphorylated tyrosine hydroxylase. Characterization of site-directed mutants at cysteine residues in 14-3-3ζ (C25, C94, and C189) by DSF and mass spectroscopy revealed covalent modification by ebselen of all cysteines through a selenylsulfide bond. C25 appeared to be the preferential site of ebselen interaction in vitro, whereas modification of C94 was the main determinant for protein destabilization. At therapeutically relevant concentrations, ebselen and ebselen oxide caused decreased 14-3-3 levels in SH-SY5Y cells, accompanied with an increased degradation, most probably by the ubiquitin-dependent proteasome pathway. Moreover, ebselen-treated zebrafish displayed decreased brain 14-3-3 content, a freezing phenotype, and reduced mobility, resembling the effects of lithium, consistent with its proposed action as a safer lithium-mimetic drug. Ebselen has recently emerged as a promising drug candidate in several medical areas, such as cancer, neuropsychiatric disorders, and infectious diseases, including coronavirus disease 2019. Its pleiotropic actions are attributed to antioxidant effects and formation of selenosulfides with critical cysteine residues in proteins. Our work indicates that a destabilization of 14-3-3 may affect the protein interaction networks of this protein family, contributing to the therapeutic potential of ebselen. SIGNIFICANCE STATEMENT: There is currently great interest in the repurposing of established drugs for new indications and therapeutic targets. This study shows that ebselen, which is a promising drug candidate against cancer, bipolar disorder, and the viral infection coronavirus disease 2019, covalently bonds to cysteine residues in 14-3-3 adaptor proteins, triggering destabilization and increased degradation in cells and intact brain tissue when used in therapeutic concentrations, potentially explaining the behavioral, anti-inflammatory, and antineoplastic effects of this drug.

Polyethylene glycol perturbs the unfolding of CRABP I: A correlation between experimental and theoretical approach.

The importance of macromolecules paves the way towards a detailed molecular level investigation as all most all cellular processes occurring at the interior of cells in the form of proteins, enzymes, and other biological molecules are significantly affected because of their crowding. Thus, exploring the role of crowding environment on the denaturation and renaturation kinetics of protein molecules is of great importance. Here, CRABP I (cellular retinoic acid binding protein I) is employed as a model protein along with different molecular weights of Polyethylene glycol (PEG) as molecular crowders. The experimental evaluations are done by accessing the protein secondary structure analysis using circular dichroism (CD) spectroscopy and unfolding kinetics using intrinsic fluorescence of CRABP I at 37 °C to mimic the in vivo crowding environment. The unfolding kinetics results indicated that both PEG 2000 and PEG 4000 act as stabilizers by retarding the unfolding kinetic rates. Both kinetic and stability outcomes presented the importance of crowding environment on stability and kinetics of CRABP I. The molecular dynamics (MD) studies revealed that thirteen PEG 2000 molecules assembled during the 500 ns simulation, which increases the stability and percentage of ß-sheet. The experimental findings are well supported by the molecular dynamics simulation results.

Choice of Force Field for Proteins Containing Structured and Intrinsically Disordered Regions.

Biomolecular force fields optimized for globular proteins fail to properly reproduce properties of intrinsically disordered proteins. In particular, parameters of the water model need to be modified to improve applicability of the force fields to both ordered and disordered proteins. Here, we compared performance of force fields recommended for intrinsically disordered proteins in molecular dynamics simulations of three proteins differing in the content of ordered and disordered regions (two proteins consisting of a well-structured domain and of a disordered region with and without a transient helical motif and one disordered protein containing a region of increased helical propensity). The obtained molecular dynamics trajectories were used to predict measurable parameters, including radii of gyration of the proteins and chemical shifts, residual dipolar couplings, paramagnetic relaxation enhancement, and NMR relaxation data of their individual residues. The predicted quantities were compared with experimental data obtained within this study or published previously. The results showed that the NMR relaxation parameters, rarely used for benchmarking, are particularly sensitive to the choice of force-field parameters, especially those defining the water model. Interestingly, the TIP3P water model, leading to an artificial structural collapse, also resulted in unrealistic relaxation properties. The TIP4P-D water model, combined with three biomolecular force-field parameters for the protein part, significantly improved reliability of the simulations. Additional analysis revealed only one particular force field capable of retaining the transient helical motif observed in NMR experiments. The benchmarking protocol used in our study, being more sensitive to imperfections than the commonly used tests, is well suited to evaluate the performance of newly developed force fields.

Structural Basis of Ca2+-Dependent Self-Processing Activity of Repeat-in-Toxin Proteins.

The posttranslational Ca2+-dependent "clip-and-link" activity of large repeat-in-toxin (RTX) proteins starts by Ca2+-dependent structural rearrangement of a highly conserved self-processing module (SPM). Subsequently, an internal aspartate-proline (Asp-Pro) peptide bond at the N-terminal end of SPM breaks, and the liberated C-terminal aspartyl residue can react with a free ε-amino group of an adjacent lysine residue to form a new isopeptide bond. Here, we report a solution structure of the calcium-loaded SPM (Ca-SPM) derived from the FrpC protein of Neisseria meningitidis The Ca-SPM structure defines a unique protein architecture and provides structural insight into the autocatalytic cleavage of the Asp-Pro peptide bond through a "twisted-amide" activation. Furthermore, in-frame deletion of the SPM domain from the ApxIVA protein of Actinobacillus pleuropneumoniae attenuated the virulence of this porcine pathogen in a pig respiratory challenge model. We hypothesize that the Ca2+-dependent clip-and-link activity represents an unconventional strategy for Gram-negative pathogens to adhere to the host target cell surface.IMPORTANCE The Ca2+-dependent clip-and-link activity of large repeat-in-toxin (RTX) proteins is an exceptional posttranslational process in which an internal domain called a self-processing module (SPM) mediates Ca2+-dependent processing of a highly specific aspartate-proline (Asp-Pro) peptide bond and covalent linkage of the released aspartyl to an adjacent lysine residue through an isopeptide bond. Here, we report the solution structures of the Ca2+-loaded SPM (Ca-SPM) defining the mechanism of the autocatalytic cleavage of the Asp414-Pro415 peptide bond of the Neisseria meningitidis FrpC exoprotein. Moreover, deletion of the SPM domain in the ApxIVA protein, the FrpC homolog of Actinobacillus pleuropneumoniae, resulted in attenuation of virulence of the bacterium in a pig infection model, indicating that the Ca2+-dependent clip-and-link activity plays a role in the virulence of Gram-negative pathogens.

Quantum Chemical Calculations of NMR Chemical Shifts in Phosphorylated Intrinsically Disordered Proteins.

Quantum mechanics (QM) calculations are applied to examine 1H, 13C, 15N, and 31P chemical shifts of two phosphorylation sites in an intrinsically disordered protein region. The QM calculations employ a combination of (1) structural ensembles generated by molecular dynamics, (2) a fragmentation technique based on the adjustable density matrix assembler, and (3) density functional methods. The combined computational approach is used to obtain chemical shifts (i) in the S19 and S40 residues of the nonphosphorylated and (ii) in the pS19 and pS40 residues of the doubly phosphorylated human tyrosine hydroxylase 1 as the system of interest. We study the effects of conformational averaging and explicit solvent sampling as well as the effects of phosphorylation on the computed chemical shifts. Good to great quantitative agreement with the experiment is achieved for all nuclei, provided that the systematic error cancellation is optimized by the choice of a suitable NMR standard. The effect of the standard reference on the computed 15N and 31P chemical shifts is demonstrated by employing three different referencing methods. Error bars associated with the statistical averaging of the computed 31P chemical shifts are larger than the difference between the 31P chemical shift of pS19 and pS40. The sequence trend of 31P shifts therefore could not be reliably reproduced. On the contrary, the calculations correctly predict the change of the 13C chemical shift for CB induced by the phosphorylation of the serine residues. The present work demonstrates that QM calculations coupled with molecular dynamics simulations and fragmentation techniques can be used as an alternative to empirical prediction tools in the structure characterization of intrinsically disordered proteins.

Structure and Functions of Microtubule Associated Proteins Tau and MAP2c: Similarities and Differences.

The stability and dynamics of cytoskeleton in brain nerve cells are regulated by microtubule associated proteins (MAPs), tau and MAP2. Both proteins are intrinsically disordered and involved in multiple molecular interactions important for normal physiology and pathology of chronic neurodegenerative diseases. Nuclear magnetic resonance and cryo-electron microscopy recently revealed propensities of MAPs to form transient local structures and long-range contacts in the free state, and conformations adopted in complexes with microtubules and filamentous actin, as well as in pathological aggregates. In this paper, we compare the longest, 441-residue brain isoform of tau (tau40), and a 467-residue isoform of MAP2, known as MAP2c. For both molecules, we present transient structural motifs revealed by conformational analysis of experimental data obtained for free soluble forms of the proteins. We show that many of the short sequence motifs that exhibit transient structural features are linked to functional properties, manifested by specific interactions. The transient structural motifs can be therefore classified as molecular recognition elements of tau40 and MAP2c. Their interactions are further regulated by post-translational modifications, in particular phosphorylation. The structure-function analysis also explains differences between biological activities of tau40 and MAP2c.

Functionally specific binding regions of microtubule-associated protein 2c exhibit distinct conformations and dynamics.

Microtubule-associated protein 2c (MAP2c) is a 49-kDa intrinsically disordered protein regulating the dynamics of microtubules in developing neurons. MAP2c differs from its sequence homologue Tau in the pattern and kinetics of phosphorylation by cAMP-dependent protein kinase (PKA). Moreover, the mechanisms through which MAP2c interacts with its binding partners and the conformational changes and dynamics associated with these interactions remain unclear. Here, we used NMR relaxation and paramagnetic relaxation enhancement techniques to determine the dynamics and long-range interactions within MAP2c. The relaxation rates revealed large differences in flexibility of individual regions of MAP2c, with the lowest flexibility observed in the known and proposed binding sites. Quantitative conformational analyses of chemical shifts, small-angle X-ray scattering (SAXS), and paramagnetic relaxation enhancement measurements disclosed that MAP2c regions interacting with important protein partners, including Fyn tyrosine kinase, plectin, and PKA, adopt specific conformations. High populations of polyproline II and α-helices were found in Fyn- and plectin-binding sites of MAP2c, respectively. The region binding the regulatory subunit of PKA consists of two helical motifs bridged by a more extended conformation. Of note, although MAP2c and Tau did not differ substantially in their conformations in regions of high sequence identity, we found that they differ significantly in long-range interactions, dynamics, and local conformation motifs in their N-terminal domains. These results highlight that the N-terminal regions of MAP2c provide important specificity to its regulatory roles and indicate a close relationship between MAP2c's biological functions and conformational behavior.

Free energy calculations on the stability of the 14-3-3ζ protein.

Mutations of cysteine are often introduced to e.g. avoid formation of non-physiological inter-molecular disulfide bridges in in-vitro experiments, or to maintain specificity in labeling experiments. Alanine or serine is typically preferred, which usually do not alter the overall protein stability, when the original cysteine was surface exposed. However, selecting the optimal mutation for cysteines in the hydrophobic core of the protein is more challenging. In this work, the stability of selected Cys mutants of 14-3-3ζ was predicted by free-energy calculations and the obtained data were compared with experimentally determined stabilities. Both the computational predictions as well as the experimental validation point at a significant destabilization of mutants C94A and C94S. This destabilization could be attributed to the formation of hydrophobic cavities and a polar solvation of a hydrophilic side chain. A L12E, M78K double mutant was further studied in terms of its reduced dimerization propensity. In contrast to naïve expectations, this double mutant did not lead to the formation of strong salt bridges, which was rationalized in terms of a preferred solvation of the ionic species. Again, experiments agreed with the calculations by confirming the monomerization of the double mutants. Overall, the simulation data is in good agreement with experiments and offers additional insight into the stability and dimerization of this important family of regulatory proteins.

Conformational dynamics are a key factor in signaling mediated by the receiver domain of a sensor histidine kinase from Arabidopsis thaliana.

Multistep phosphorelay (MSP) cascades mediate responses to a wide spectrum of stimuli, including plant hormonal signaling, but several aspects of MSP await elucidation. Here, we provide first insight into the key step of MSP-mediated phosphotransfer in a eukaryotic system, the phosphorylation of the receiver domain of the histidine kinase CYTOKININ-INDEPENDENT 1 (CKI1RD) from Arabidopsis thaliana We observed that the crystal structures of free, Mg2+-bound, and beryllofluoridated CKI1RD (a stable analogue of the labile phosphorylated form) were identical and similar to the active state of receiver domains of bacterial response regulators. However, the three CKI1RD variants exhibited different conformational dynamics in solution. NMR studies revealed that Mg2+ binding and beryllofluoridation alter the conformational equilibrium of the ß3-α3 loop close to the phosphorylation site. Mutations that perturbed the conformational behavior of the ß3-α3 loop while keeping the active-site aspartate intact resulted in suppression of CKI1 function. Mechanistically, homology modeling indicated that the ß3-α3 loop directly interacts with the ATP-binding site of the CKI1 histidine kinase domain. The functional relevance of the conformational dynamics observed in the ß3-α3 loop of CKI1RD was supported by a comparison with another A. thaliana histidine kinase, ETR1. In contrast to the highly dynamic ß3-α3 loop of CKI1RD, the corresponding loop of the ETR1 receiver domain (ETR1RD) exhibited little conformational exchange and adopted a different orientation in crystals. Biochemical data indicated that ETR1RD is involved in phosphorylation-independent signaling, implying a direct link between conformational behavior and the ability of eukaryotic receiver domains to participate in MSP.

Exploring the binding pathways of the 14-3-3ζ protein: Structural and free-energy profiles revealed by Hamiltonian replica exchange molecular dynamics with distancefield distance restraints.

The 14-3-3 protein family performs regulatory functions in eukaryotic organisms by binding to a large number of phosphorylated protein partners. Whilst the binding mode of the phosphopeptides within the primary 14-3-3 binding site is well established based on the crystal structures of their complexes, little is known about the binding process itself. We present a computational study of the process by which phosphopeptides bind to the 14-3-3ζ protein. Applying a novel scheme combining Hamiltonian replica exchange molecular dynamics and distancefield restraints allowed us to map and compare the most likely phosphopeptide-binding pathways to the 14-3-3ζ protein. The most important structural changes to the protein and peptides involved in the binding process were identified. In order to bind phosphopeptides to the primary interaction site, the 14-3-3ζ adopted a newly found wide-opened conformation. Based on our findings we additionally propose a secondary interaction site on the inner surface of the 14-3-3ζ dimer, and a direct interference on the binding process by the flexible C-terminal tail. A minimalistic model was designed to allow for the efficient calculation of absolute binding affinities. Binding affinities calculated from the potential of mean force along the binding pathway are in line with the available experimental estimates for two of the studied systems.